o 
o 

(N 



o 



T3 
O 

o 



> 
in 

m 
o 

o 

-I— > 



i 

-a 
c 

o 
o 



X 
S3 



Non-linear screening of spherical and cylindrical colloids: the case of 1:2 and 2:1 

electrolytes 
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From a multiple scale analysis, we find an analytic solution of spherical and cylindrical Poisson- 
Boltzmann theory for both a 1:2 (monovalent co-ions, divalent counter-ions) and a 2:1 (reversed 
situation) electrolyte. Our approach consists in an expansion in powers of rescaled curvature l/(«a), 
where a is the colloidal radius and 1/k the Debye length of the electrolytic solution. A systematic 
comparison with the full numerical solution of the problem shows that for cylinders and spheres, 
our results are accurate as soon as na > 1. We also report an unusual overshooting effect where the 
colloidal effective charge is larger than the bare one. 
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I. INTRODUCTION 

Almost a century ago, the work of Gouy pj followed 
by that of Chapman |2( , has established the foundations 
of the mean-field treatment of the electric-double layer 
(Poisson-Boltzmann theory). This approach served as a 
basis for computing microionic correlations in a homo- 
geneous electrolyte and later led to the prominent 
DLVO theory of colloidal interactions Q. An essential 
notion in this context is that of charge renormalization 
IE IE IE 13 : a t large distances, the electrostatic sig- 
nature of a charged body (with charge Zb are ) in an elec- 
trolyte takes the same form as that of an effective macro- 
ion with a suitable effective charge Z c g, the latter object 
being treated within linearized Poisson-Boltzmann the- 
ory. Only for small Zb ar c do effective and bare parame- 
ters coincide (weak coupling limit). In general, one has 
\Z c g\ <C |^bare| which reflects the non- linear screening 
effect of the electric double layer around a colloid [Io| . 
This non-linear regime, beyond the weak coupling limit 
but below the couplings that would invalidate the mean- 
field assumption underlying the approach, is precisely 
that which is relevant for colloids (see e.g. the discussion 
in references 0,E1). 

Recently, analytical expressions have been obtained, 
within Poisson-Boltzmann theory, for the effective 
charges of spherical and cylindrical macro-ions |12| . 
These predictions for a unique macro-ion immersed in 
an infinite sea of monovalent electrolyte with inverse De- 
bye length k, are exact up to (na)^ 1 corrections, where 
a is the radius of the macro-ion. For practical purposes, 
the predictions are accurate as soon as kcl > 1. In this 
paper, we consider the situation of spherical and cylin- 
drical macro-ions in a charge asymmetric electrolyte with 



both monovalent and divalent micro-ions. The asym- 
metry of electrolyte has noticeable consequences on the 
structure of the electric double-layer and the case of 2:1 
electrolytes (i.e. with divalent co-ions and monovalent 
counter- ions) turns out to differ much from the 1:2 sit- 
uation (monovalent co-ion, divalent counter-ion). Our 
analytical results -obtained from a multiple scale tech- 
nique neglect 0(ko)~ 2 corrections for the electro- 
static potential and conversely, 0{k,cl)~ 1 terms for effec- 
tive charges. By an explicit comparison with the numer- 
ical mean-field results, they will be shown to be precise 
whenever na > 1, as was the case in I n section 
ITU the general method will be presented, and the elec- 
trostatic potential obtained. The results concerning ef- 
fective quantities will be given in sections IIIII and IIVI 
Conclusions will be drawned in section Ivl 



II. QUASI-PLANAR SOLUTION TO 
POISSON-BOLTZMANN EQUATION FOR 2:1 OR 
1:2 ELECTROLYTES 

A. 2:1 electrolyte 

We consider a cylindrical (j — 1) or spherical (j = 2) 
colloid of radius a with surface charge density ea > im- 
mersed in an electrolyte with co-ions (resp. counterions) 
of valency z\ (resp. Z2) and numeric density n\ (resp. 712). 
Let us analyze in some detail the case z\ = 2, zi = — 1, 
hereafter referred to as 2:1. 



As 



usual, we define the Debye length k 

,2\ 
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(A-kIb J2i nizf)^ 1 / 2 = (127rrt 2 ^s) 1 ^ 2 , the reduced electri- 
cal potential y — 0eip and a* = ^irlscm. Here, Ib de- 
notes the Bjerrum length, defined from the permittivity 
X of the suspending medium and the inverse temperature 
(3 = l/(kT) as Ib = (3e 2 /x- Using the method of multi- 
ple scales, closely following ^Eli the Poisson-Boltzmann 



2 



equation 

1 d 

ri dr 

can be cast into 



dr 



= -ATrl B n 2 {e- 2y ~e v ) (2.1) 



d 2 y 2ed 2 y ej dy 2 ^ 2 V ^ ®V 
dx\ dx\dx2 xi dx\ dx\ xi dxi 

1 



-(e-^-ey) (2.2) 



with boundary conditions 
dy dy 



edxi 0x2 



lim 



X\ — >OC,X*2 — >00 



X\ — 0,X2 — 1 

dy dy 
edxi dx2 



(2.3a) 
(2.3b) 



Here, we have defined e = (na)^ 1 , x± = n(r — a) and 
x% = r/a. We seek for a solution as an expansion in 
powers of e which is supposed to be a small parameter: 

V = Vo + m H . 

The equation for the zeroth order term is Poisson- 
Boltzmann equation for a planar interface 



d 2 y 
dx\ 



1 



(e 



-2yo 







(2.4) 



which has been solved by Gouy in his pioneering work [lj 
(see also Grahame 14]). The solution reads 



2/o(zi,Z2) = In 1 + 



6q 



{1-qf 



(2.5) 



with the short-hand notation q = t(x2)e~ Xl , which will 
be used extensively in the following. Here t(x2) is a func- 
tion of X2 which appears as a constant of integration 
(with respect to x%) since in Eq. H2.4f> the variable X2 
does not appear. As explained in [13} this function is de- 
termined by the requirement that the non-homogeneous 
part of the differential equation for the next order, yi, 
decays faster than e~ Xl when x\ — > oo. The equation for 
yi reads 



d 2 Vl 
dx\ 



1 



(2e- 2yo + e V0 ) yi 



2d 2 y Q 
dx\dx2 



3 dy 
X2 dx\ 



(2.6) 



The requirement that the r.h.s. of Eq. (|2.6(l decays faster 

than e~ Xl leads to t(x2) = Ax^^ 2 with A a constant of 
integration. We therefore have 



Ax. 



-J/2 -2 



(2.7) 



Notice that the situation is exactly the same as in the 
1:1 electrolyte case the zero order solution in the 
quasi-planar approximation is obtained from the planar 
solution with the replacement of the constant of inte- 
gration A by Ax 2 ^ 2 . Actually, this is a general result 



for any type of electrolyte since the r.h.s. of Eq. I|2.6() 
does not depend on the microscopic constitution of the 
electrolyte and when x\ — > oo for any electrolyte the be- 
havior of j/o will be given by the Debye-Hiickel solution: 
est x t{x?) exp(— x\). 

The constant of integration A can be expressed as a 
function of the surface charge density a* by enforcing 
the boundary condition l|2.3a|l at the dominant order 



dyo 
dxx 



(2. 



where we have put s = ea* . This gives a third order 
equation for A 



6A(l+A) 



(1- A){A 2 +AA+1) 



(2.9) 



Its physical solution (which vanishes when s — * 0) can be 
written as 



A = l 

s 



-2- S + 2 3 / 2 (2 + S + S 2 ) 1 / 2 cos(£) 



with 



cos 



-4 - 3s - 3s 2 - s 3 
V2(2 + s + s 2 ) 3 / 2 



(2.10) 



(2.11) 



This constant has also been computed in the study of 
the plan ar interface effective charge (6A — 47rs c ff) in 
Ref. [l5| , although it is presented there in a slightly (but 
completely equivalent) form. 

Replacing the explicit expression l|2.5|l for yo into 
Eq. (|2.t)|) gives for the order-one term y\ the following 
equation 



yi 



d 2 y x 1 + 27g 2 + 16g 3 + 27q 4 + q 6 
~dxj (l-q) 2 (l+4q + q 2 ) 2 

12j q 2 (q 3 + 3q 2 + 3q - 1 



x 2 (l-9) 2 (l + 4<7 + g 2 ) 5 



(2.12) 



Using the variable q instead of x\ and performing the 
change of function yi(xi,x 2 ) = f{q)/[0--q){^ + ^q + q 2 )] 
yields a second order linear differential non-homogeneous 
equation for f(q) with polynomial coefficients in q whose 
associated linear homogeneous equation has the simple 
solution f(q) — q(q + 1), therefore allowing to find the 
complete solution to the non-homogeneous equation us- 
ing the traditional method of "variation of the constant" . 
After some tedious but otherwise straightforward calcu- 
lations, we find the solution satisfying the appropriate 
boundary condition (|2.3b|) at infinity, 



yi{xi,x 2 ) 



k(x 2 )q(q+l)- ^~q 2 (q 2 +9q 



(l-q)(l + 4q + q 2 ) 

(2.13) 

Again there is a function k(x2) that appears as a con- 
stant of integration with respect to x\ since there are 



3 



no derivatives of y\ with respect to x-i in Eq. (|2.6|l . This 
function hix^) is determined [l3| by the requirement that 
the non-homogeneous part of the equation for the next 
order term yi decreases faster that e~ Xl when x\ — > oo 



d 2 m 

dx\ 



3 dyo 2 _^m_ + J_^yi 

X2 dx2 dx\dx2 X2 dxi 
+y 2 1 {e yo -4e- 2yo ) 



&i£ (2-14) 



= o(e- ai ) 



This gives k{x2) = c\ + 3j(j — 2)/(4x2) with c% another 



constant of integration, so finally the order-one solution 
is 



yi(xi,x 2 ) = 



Cl 



3j(j-2) 
4x 2 



(l-q)(l+4q + q 2 ) 



(2.15) 

Applying the boundary condition i|2.3a|) to the next order 
in e gives the equation d Xl yi + d X2 y \ Xl=0 X2=1 = and 
subsequently determines the constant of integration 



Cl 



2A 6 + 12A 5 + A 4 {34 + 3j) + 2A 3 (-88 + 3j) + 6A 2 (-7 + 3j) + 2A(34 + 3j) + 3(2 + j) 
1 4(l + 2^ + 6A 2 + 2A 3 + ^l 4 ) 



(2.16) 



r 



The quantity A is given by Eq. (|2.1UI) . Both constants A 
and Ci are related to the effective charge of the colloid and 
therefore carry important physical information about the 
system. Let us notice that at saturation er — > oo, they 
take simple values: yl sat = 1 and c{ at = — 8)/4. 



B. 1:2 electrolyte 

The quasi-planar approximate solution of the Poisson- 
Boltzmann equation for the case Z\ = 1 and z 2 = —2 (1:2 
electrolyte) follows from similar calculations. We only 
report the results. The zero order term yo reads 



Vo(xi,x 2 ) 



In 1 



(q+l) 2 

with q given by Eq. (|2.7|l and the order-one term is 



(2.17) 



yi{xi,x 2 ) = 



-9(9-1) 



ci 



3j(j-2) 

4X2 



+ 2^1 {q -9<? 



(l+q)(l-4q + q3) 



Note that the solution for the 1:2 case is simply obtained 
from the one for the 2:1 case by a global change of sign 
and by replacing q by —q. 

The constant of integration A is again a solution 
of a third order equation which can be obtained from 
Eq. (|2.9|l by a global change of sign and by replacing A 
by —A. However, the physical solution is not the same 
as in the 2:1 case, and now takes the form 



.4 



2 + a + 2 3 / 2 (2 - s + s 2 ) 1 ' 2 cos (^p^ 
with 9 given by 



(2.19) 



-4 + 3s - 3s 2 + s 6 



V2(2- 



2)3/2 



(2.20) 



(2.18) The constant of integration for the order-one term is here 



Cl 



2A 6 - 12A 5 + A 4 (34 + 3j) - 2A 3 (-88 + 3j) + 6A 2 (-7 + 3j) - 2A(34 + 3j) + 3(2 + j) 
' 4(1 - 2A + 6A 2 - 2A 3 + A*) : 

I 



(2.21) 



with A given by Eq. I|2.19|l . The saturation (s — > oo) 
values of these constants are now different. We have 
A sat = 2 - V3 and cf* = - j(28 + 3j - 24V3)/4. 



Comparison between analytical and numerical 
potential profiles 



Gathering results, we obtain up to corrections of order 
1/(ko) 2 , y(r) = yo{f) + ( K a) _1 2/i( r )j where the auxiliary 
functions yo and y\ are given by Eqs. I|2.5|) and l|2.15ll in 
the 2:1 case, and by Eqs. (|2~T7|) and (|2~T%|) for 1:2 elec- 



4 




FIG. 1: Reduced electrostatic potential y(r) as a function 
of rescaled distance for a spherical macro-ion in a 1:2 elec- 
trolyte. The continuous curve shows the numerical solution 
of the problem and the crosses indicate the values found from 
Eq. Ij2.17ft and l|2.18ft . The inset shows the same data on a 
linear-log scale. Here, ko = 2 and the reduced bare charge is 
very high: Z b3bIB l B /a = 2000. 



trolytes. It is instructive to compare the resulting pre- 
dictions to the numerical solution of Poisson-Boltzmann 
theory, obtained following the method of Ref. ^5] . Fig- 
ures n and |21 show that already for na = 2, the agree- 
ment is good. Although the potential at contact y{a) 
is predicted accurately, we observe that our theoretical 
expression slightly underestimates the potential. A simi- 
lar trend will be observed for effective charges -again for 
spheres- in section ITTT1 In cylindrical geometry, a slight 
overestimation may be found in the 1:2 case. 

The parameters in figures ^ and |21 are chosen to be in 
the non- linear saturation regime Zbarc 3> a/^s- It is in- 
teresting to notice that the relative error of our analytic 
solution from the numerical one in the cases presented in 



figures n and |21 is of order 3% that is of order (a 



710. 



We have also studied the linear regime, Zbare sma ll, an( i 
in this case the error is larger of order 25%, i.e. of or- 
der (ko)~ 2 (remember that in our analytical solution we 
neglect terms of order (ko)~ 2 ). We have also computed 
the relative error for other values of na and the trend is 
general: in the linear regime the relative error is of or- 
der (na)" 2 but for the non- linear saturation regime the 
situation improves and the error is reduced by a factor 
10. This makes our analytic solution practical since ex- 
perimental situations are often in the saturation regime 
where our solution is more accurate. 



III. EFFECTIVE CHARGES 



A. Spheres 

The far field r — > oo behavior of the solution y(r) = 
yo(r) + eyi(r) + 0(e 2 ), obtained in the last section is 

y(r) ~ Ae-< r - a U-\' 2 + +0(e 2 ). (3.1) 



FIG. 2: Same as Fig. 0in a 2:1 electrolyte. Here, na = 2 and 
the reduced bare charge is: Zbarc is/a = 34. 



With this expression, we can deduce the effective charge. 
For a spherical macro-ion (j = 2) of radius a and charge 
Z e g , the solution of linearized Poisson-Boltzmann theory 
(also referred to as Debye-Huckel theory) V 2 y = n 2 y 
reads 



y(r) = 



Z eB l B e -^-") 



na 



(3.2) 



By comparison with expression (|3.1(l we conclude that 
the effective charge is given by 



A 



6na + 6 + c\ + O 



(3.3) 



The coefficients A and ci are given by Eqs. 12.10J1 and 
$FWi (2:1 electrolyte) or Eqs. fFTty and lt!T2T)l (1:2 elec- 
trolyte) in terms of the bare charge .Zbare by reporting 
s = eZhaiJ-B I 'a,. Figures 01 and 0] compare the above an- 
alytical predictions to the effective charge obtained from 
the far field behavior of the numerical solution of Poisson- 
Boltzmann theory, obtained as explained in 16]. The 
agreement is satisfying, and improves upon increasing 
na, as was anticipated. 

One may readily check from Eq. (|3.3|l that in the limit 
■Zbarc — ► 0, -Zcff/Zbarc — > 1. Effective and bare param- 
eters coincide in the weak coupling limit, as it should 
(see the dashed lines in Figures |31 and • In the other 
limit where Zbarc — > oo, we observe the saturation picture 
common to several mean-field theories 0, 0] : the effec- 
tive charge goes to a plateau value, that only depends 
on two dimcnsionless quantities, a/ls and na. The effec- 
tive charge at saturation (obtained when s — > oo) takes 
a simple expression. For a 2:1 electrolyte 



7 sat Ib 
a 



= 6na + 7 + 



(3.4) 



and for a 1:2 electrolyte 



7 sat 1 b 

a 



2-Vs) kna - 11 + 12\/3 + O (e)](3.5) 



1.608 na 



2.623 + ( — 

na 



5 



10 




r ' 1 ■ 1 ■ 1 — 1 

20 40 60 

Z b„re l B /Cl 

FIG. 3: Effective vs bare charge for a spherical macro-ion in 
a 1:2 electrolyte (i.e. monovalent co- ions/divalent counter- 
ions). The open circles are obtained from the full non-linear 
Poisson Boltzmann theory, while the continuous curve corre- 
sponds to the analytical prediction given by Eqs. 13.31 . The 
dashed line has slope 1 and shows the initial linear regime 
for weak charges. The salinity conditions are here such that 
na = 3 where a is the sphere radius. 
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FIG. 4: Same as Figure |3 for a 2:1 electrolyte (divalent co- 
ions, monovalent counter-ions). As indicated, the main graph 
corresponds to na — 2 while the inset shows results for na — 5. 
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FIG. 5: Effective charge at saturation in a 1:2 electrolyte for a 
spherical colloid. The line shows the prediction of Eq. 13.51 . 
The circles again correspond to the numerical resolution of 
Poisson-Boltzmann theory, and the dotted line between them 
is a guide to the eye. The inset shows the same data in a 
log-linear scale. 

Z s ^l B /a = 4ko + 6 holds for m > 1 [T^|. 



40 




FIG. 6: Same as Figure|^]for a 2:1 electrolyte. The line shows 
the prediction of Eq. 13.41 . 



We note that the conditions of Fig. ^ are those of satu- 
ration. 

These expressions are tested against the numerical 
data in Figures |S] and The agreement is good for 
na > 1. In Figs. [S] and El an inset has been added to 
show the regime of low na values where our approach 
breaks. In this limit, the observed divergence of 
means that the bare charge is no longer renormalized. 
As happens for monovalent electrolytes [5| , the saturated 
effective charge is a non-monotonous function of na, that 
reaches its minimum for na ~ 0.3. In the latter 1:1 case, 
we recall for completeness that the asymptotic expansion 



B. Cylinders 

For an infinite long cylindrical colloid, j = 1, with lin- 
ear charge density eA, the far-field solution l|3.1(l should 
be compared to the one obtained from Debye-Hiickcl the- 
ory, 



y(r) 



2l B \ cS K (nr) 
na K\(tia) 



(3.6) 
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-«(r— a) 
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FIG. 7: Effective line charge density as a function of its bare 
counterpart, for a rod-like macro-ion. The line shows the 
prediction of Eq. 13. 71 and the symbols stand for the "exact" 
numerical values. The main graph and the inset correspond to 
the same salinity conditions na — 3, where a is the cylinder's 
radius. 



FIG. 8: Saturated effective line charge of an infinite cylinder 
with radius a in a 1:2 electrolyte. Line: Eq. 13.91 and sym- 
bols: numerical solution. The inset is a log-linear plot. The 
arrow indicates the value \/3/(27r) ~ 0.275, obtained in the 

kb -> o limit mm. 



where Kq and K\ are modified Bessel functions. We 
conclude that the effective line charge density is given by 



3na 



£i 
2 



O 



(3.7) 



The explicit expression of the effective linear charge den- 
sity in terms of the bare linear charge density Abarc is 
obtained by reporting s = 2els Abarc in the analytical ex- 
pressions for A and c\ given in Eqs. p.lOfl and (|2.16() (2:1 
electrolyte) or Eqs. l(2~TT?| and (|2~2"Tf (1:2 electrolyte). 
Figured shows the accuracy of our analytical expression, 
that turns out to be slightly better in the 2:1 situation 
than in the 1:2 case (the reverse observation follows from 
inspecting Figures and EJ) ■ 

The effective charges at saturation are, for a 2:1 elec- 
trolyte, 



7 

4 



and for a 1:2 electrolyte, 



fflB 



(2-V3) 



3ko 







11 

T 



3^3 + O(e) 



(3.8) 



(3.9) 



0.804 na + 0.655 + 



These simple expressions are plotted in Figures [H] and [51 
together with their counterparts obtained from the so- 
lution of Poisson-Boltzmann theory, shown by symbols. 
When converted into effective surface charge densities 
dig, the previous results yield, up to (na)^ 1 corrections 

4iral B <7 s c s = 6na + 7 (2 : 1, spheres) (3.10) 
7 

= 6na + - (2 : 1, rods), (3.11) 



while in the planar case, one gets Anal b (fig = 6/ta + 
0. The increase of the zeroth order term (0, 7/2, 7) as 
the dimensionality of the object increases, reflects the 
concomitant weaker range of the bare Coulomb potential 
(— r in ID, — logr in 2D, 1/r in 3D. . . ). Indeed, a weaker 
Coulomb contribution leads to a weaker screening, hence 
a higher effective charge. A similar argument therefore 
explains the increase of effective charges with na. For 
completeness, we also give the 1:2 results 



*B <7 cfT 



(2-V3) 



6/ia- 11 + 12\/3 



6na — + 6V3 



(spheres) (3.12) 
(rods), (3.13) 



and the same argument as above equally applies here. 

It may be observed in Figures |S] and [5] that in cylin- 
drical geometry, the saturated effective charge does not 
diverge in the limit na — > 0, as was the case for spheres 
due to entropic reasons (Boltzmann "beats" Coulomb in 
this limit). With infinite cylinders, this is no longer the 
case (a two dimensional situation is more favorable to 
Coulomb than a 3D one), and A|g reaches a constant 
value for small na (see the inset of Figures |S] and |5J) . 



C. An overshooting effect 

Although the effect is not very marked, it may be ob- 
served in Fig. 01 that Z e s as a function ^bare has an 
inflexion point, so that Z c r > Zbare in a given charge 
range (0 < Z c r < lOa/Zs, where not only our predic- 
tion but also the symbols showing numerical data lie 
above the dashed line). This was unexpected since with 
a monovalent (1:1) electrolyte, the effective charge is al- 
ways smaller than the bare one. This overshooting effect 



7 



c< 10 




FIG. 9: Same as Figure for a 2:1 electrolyte. The line 
corresponds to Eq. 13.81 



occurs for spherical but also for rod-like macroions. It 
requires a 2:1 salt for which the divalent co-ions are ex- 
pelled from the vicinity of the macro-ion, which leads 
to a much weaker screening than in the reverse 1:2 sit- 
uation. That this effect is able to impose Z c g > Zb aic 
is however surprising, and in order to check its robust- 
ness, we also investigated numerically more asymmetric 
electrolytes. Figure ITU1 shows that for a 5:1 salt, the over- 
shooting is more pronounced and that for Z e g ~ lOa/ls, 
the effective charge may be twice bigger than the bare 
one. 

Additional insight into this unexpected overshooting 
effect may be obtained from our analytical expressions 
for the effective charges Eqs. <|3.3[) and (|3.7|) . For a small 
surface charge density s, in the 2:1 case, the quantities A 
and Ci involved in the expressions of the effective charges 
have a Taylor expansion of the form 



A 

ci 
ci 



s s s 
6 + 18 ~ 216 
9 7s 13s 2 

~4 ~ ¥ 
14s 



+ 0(s 4 ) 
10s 3 



-6 



24 
13s 2 



9 

20s 3 



9 



(3.14a) 

fO(s 4 ); rods(3.14b) 

+ 0(s ); spheres 

(3.14c) 



For rods, this gives the following behavior of the effective 
linear charge density for small bare charge 



AnffZ 



cff'S 



Abarc^B + " 



, na+ — ] +0(s 4 ) 

72 v i2 ; V ' 



V 

6 V " 6 y 

(3.15) 

In this case s = 2Ab ar e^B(«a) 1 . For spherical colloids, 
the behavior for small bare charges is 



7 Ib 7 IbZs 2 ( 7\ s 3 
Z eS -=Z hale -+— Ua-- -- 



f)+<V> 

(3.16) 



with now s = Z^ aie lB/(Ka 2 ). Remembering that our an- 
alytical solution is valid for large values of na, we notice 
that in both cases the coefficient of the term of order two 
in the bare surface charge density s is positive. This im- 
plies that the effective charge will become larger than the 
bare charge in a certain intermediate regime of values of 
the bare charge when non-linear effects start to become 
important (being nevertheless far form the strongly non- 
linear saturated regime where the effective charge satu- 
rates). Let us mention that this overshooting effect is 
also expected for a planar geometry, since in that case 
the effective charge is essentially given by A. 

In contrast with this, in the case of a 1:2 electrolyte the 
Taylor expansions of A and c\ are similar to Eqs. (|3.14|) 
with the formal replacement s — > — s and A — > —A, in 
particular the sign of the order s 2 in A changes sign. This 
gives the following behavior for the effectives charges 
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for spheres. The coefficient of s 2 has changed sign with 
respect to the 2:1 case. This coefficient is now nega- 
tive and this implies that the effective charge will remain 
smaller that the bare charge, thus no overshooting effect 
for the 1:2 electrolyte case. 

It is interesting to mention that for a symmetric 1:1 
electrolyte, the small bare charge behavior reads 01 



Acff^S — Abare^fl f 



- +0(x>) (3.19) 



for rods with x r = Abaro^s/I^a + (1/2)] and 

r 



ry g ry & S 

^eff — ^bare ^~ 

a a 2 



+ 0(x>) (3.20) 



for spheres with x s = Z\, al J,B / '[2a(na +1)]- Interest- 
ingly, there is no term of order two in the bare charge as 
opposed to the asymmetric electrolytes cases. The first 
correction to the linear term is of order three and it is 
negative. The effective charge will be smaller than the 
bare charge: no overshooting effect here either. 

The first corrections to the linear theory are of order 
two in the bare charge for asymmetric electrolytes and 
with the sign of z^n a , z a being the valency of species 
a and n a its density. On the other hand, for symmetric 
electrolytes, the first correction is of order three in the 
bare charge. This important difference between symmet- 
ric and asymmetric electrolytes also appears in others 
contexts, namely in the study of the contributions due to 
correlations to the effective charge, in a framework going 
beyond the mean field approximation [rH [l8L Il9| . 

Finally, we mention that recent HNC integral equa- 
tion computations for the same systems as investigated 
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FIG. 10: Illustration of the overshooting effect, for a spherical 
colloid in a 5:1 electrolyte, with na = 2. The dashed line has 
slope 1 and the inset (a zoom on the bottom left corner) shows 
that for small charges, Z e g is a convex-up function of Zbare- 



FIG. 11: Effective potentials at saturation following from Eqs. 
14. Ill and 14. 21 . as a function of na. The limiting values for 
na — > oo are shown by the dotted lines. 



here confirm the validity of the overshooting effect |20| . 
Explicit comparisons with our predictions are under way 



IV. COLLOIDS AS CONSTANT POTENTIAL 
OBJECTS 

Colloids are usually highly charged so that their effec- 
tive charge -within mean-field- is saturated and therefore 
independent of the bare one. Yet, the bare charge is often 
not large enough to meet the region of high micro-ionic 
electrostatic correlations where the mean field approach 
would break down 0, El ■ This remark has led to the 
proposal to consider highly charged colloids as objects of 
fixed effective potential in the case of a 1:1 electrolyte 
Similar considerations may be put forward here. 
From the analysis of section IIIII the surface potentials 
y = eip/(kT) associated with effective charges read, for 
spheres 



For rod-like polyions, we get 



ye!* = 6 



i 



1 + na 



(2:1) 



(4.1) 



& = 6(2 - V3) + (2 - V3) 12 /f 17 ; (1 : 2) 

1 + na 

The important point is that the Ka dependence is very 
weak for Ka > 1, which reinforces the picture of constant 
potential objects. One may therefore consider a highly 
charged sphere as an effective body of potential 6kT/e 
or 6(2 — y/S)kT/e depending on 2:1 or 1:2 asymmetry, 
irrespective of physico-chemical parameters. In a 1:1 salt, 
one gets a value 4fcT/e jg. It is natural to find this 
quantity in between the two bounds 6kT/e and 6(2 — 
\/3)kT/e, since screening is all the more efficient as the 
valency of counter-ions is large and the valency of co-ions 
is low (in absolute values). 



sat 
VcS 



7 \ KnUa) 

6+^- hrH; (2:1) 
ZKa J Ki(Ka) 



(4- 



& - 6(2-^ ) + (2-V3) 6 ^- 11/2 ^ a i ; (1: 

Ka Ai(rea) 

Expressions (|4.1|l and (|4.2|) are plotted in Fig. ^2 

V. CONCLUSION 

In conclusion, we have found an analytical solution 
of cylindrical and spherical Poisson-Boltzmann equation 
in asymmetric 1:2 and 2:1 electrolytes. Our approach 
amounts to performing a curvature expansion, and ne- 
glects corrections of order l/(rea) 2 for the electrostatic 
potential. For Ka > 1, the corresponding solution and as- 
sociated effective charge are in excellent agreement with 
their counterparts obtained from the full numerical reso- 
lution of the problem. 

Our multiple scale analysis relies on the possibility 
to solve analytically the planar problem (correspond- 
ing to Ka — > oo). Since for a n:m electrolyte (where n 
and m respectively stand for the valency of co-ions and 
counter- ions), this solution is only known explicitly for 
n/m — 1, 1/2 and 2, we focused on 1:2 and 2:1 salts. The 
monovalent 1:1 situation has been investigated in |T^ . ll3| . 
For a given salinity k, non-linear screening is more effi- 
cient with divalent than with monovalent counter-ions. 
Accordingly, we always found higher effective charges in 
the 2:1 than in the 1:2 situation. Surprisingly, we found 
that 2:1 screening is even able to drive the effective charge 
in a regime where it is higher than the bare one. This 
overshooting effect happens in an intermediate charge 
range, since when Zb ar c is large enough, Z c g saturates. 
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